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Abstract 

P 

The information entropy of a nuclear density distribution is calculated for 
a number of nuclei. Various phenomenological models for the density distribu- 
tion using different geometry are employed. Nuclear densities calculated within 
various microscopic mean field approaches are also employed. It turns out that 
the entropy increases on going from crude phenomenological models to more so- 
phisticated (microscopic) ones. It is concluded that the larger the information 
entropy, the better the quality of the nuclear density distribution. An alternative 
approach is also examined: the net information content i.e. the sum of infor- 
mation entropies in position and momentum space S r + Sk- It is indicated that 
S r + Sk is a maximum, when the best fit to experimental data of the density and 
momentum distributions is attained. 



1 Introduction 

Information-theoretic methods have started recently to be important for the study of 
quantum mechanical systems, [I], fj, |3], |, ^|, || ^, EL mainly through the application of 



1 



the maximum entropy principle (MEP) |T0| , |TTJ. This is done by employing a suitably 
defined quantum entropy that measures the lack of information associated with the 
distribution of a quantum state over a given known basis. The MEP provides the least 
biased description consistent with the available relevant information. 

In ]12] Ghosh, Berkowitz and Parr (GBP) defined, within the ground-state density 
functional framework, several local quantities of thermodynamic-like character in order 
to characterize a non-homogenous many-electron system. Specifically, they introduced 
the concept of a local temperature corresponding to the electronic motion and an 
information entropy S associated with the electron distribution of an atomic system. 
An interesting result |13| is that S increases with the quality of the wavefunction. 

Gadre and Bendale (GB) || U investigated a different MEP procedure for atomic 
systems. They studied basis-set effects on the entropy sum S r + Sk, where <S r is the 
position-space entropy and Sk is the momentum-space entropy. They found that this 
sum increases with the basis-set quality, so that S r + Sk can be used as a measure of 
the quality of the electronic density distribution. 

A recent attempt to apply information theory to nuclear systems is ref. |§ where 
the sum S r + Sk has been calculated employing two phenomenological models for the 
nucleus. An interesting result was that the same functional form S = aN + bNlnN 
for the entropy as function of the particle number N (electrons or nucleons) holds 
approximately for atomic and nuclear systems. 

The nucleon distributions in nuclei is of fundamental importance in our under- 
standing of nuclear properties. There are several functional forms which are used for 
the description of the nuclear distributions. Their parameters are usually fitted to 
the available empirical data for the mean square radii. Moreover, density dependent 
Hartree-Fock (DDHF) |14| calculations with Skyrme or Gogny type interactions or the 
relativistic mean field (RMF) Jl5[ calculations yield fairly reliable ground state nuclear 
properties including the densities. 

In the present paper we attempt to apply the approaches of GBP and GB as a 
criterion of the quality of nuclear density distributions derived according to various 
nuclear models. We hope that this will lead to a method to choose the best nuclear 
density in different situations of practical interest, although as the models increase in 
sophistication, it turns out that the value of S stops to increase and tends to saturation. 

The paper is organized as follows: In section 2 we present the formalism of the 
information entropy according to GBP . In section 3 the numerical results are shown 
and discussed. In section 4 an alternative approach (GB) i.e. the sum S r -\- iS/j is 
discussed. In section 5 numerical results for GB are presented and finally section 6 
summarises our conclusions. 
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2 The formalism of information entropy in phase- 
space 



In [|T|, [L3| Ghosh, Berkowitz and Parr introduced a phase-space distribution f(f,p) 
associated with the ground-state density p(f) of a N-electron system. The distribution 
f(f,p) satisfies the relations 



p(r) = J dpf(r,p) 
dfp(f) = N (1) 



t(r;p) = ~ / dpp 2 f{r,p) 

where £(f ; p) is the kinetic energy density. 

Then, they defined an entropy density s(f) and entropy S associated with the 
electron density p(f) in terms of f(f,p) 

8 (f) = -kf dpf(lnf - 1) (2) 



S = J dfs(f) (3) 



and obtained the most probable distribution f(r,p) by maximizing this functional 
subject to the constraints of correct density p(f) and correct kinetic energy density 
t(f; p). 

The phase-space distribution function f(f,p) is given by a local Maxwell-Boltzmann 
distribution law: 

f(f, p) = ((3(f) /2n) 3 / 2 p(f)exp(~p 2 (3(f)) (4) 

where j3(f) = l/kT(f). The local temperature T(r) is defined by the ideal-gas expres- 
sion for the kinetic energy: 

t(f;p) = (3/2)p(f)kT(f) (5) 
Substituting eq. (4) into eq. (3) yields the Sackur- Tetrode equation [|12| for the entropy: 

S = kjp[^-lnp + -ln(2nkT)}dr (6) 

The entropy density can also be rewritten in the form: 

s(r) = -kp[ln(t/t TF ) + c] (7) 
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where txF — CkP 5 ^ 3 is the Thomas- Fermi kinetic energy density. For convenience, r 
and p are dropped from the argument of t. The constants c and c k are given by: 

c=^ + ln(4nc k /3), c k = ^(3vr 2 ) 2 / 3 (8) 

In it was found that the entropy increases with the quality of the electron 
wavefunction. As stated in |12[, one can omit all considerations for f(f,p) and simply 
postulate the existence of (7) and (3). Thus we postulate these expressions for another 
many fermion system i.e. the nucleus. In the present article, we attempt to answer the 
question: Does this S increase with the quality of nuclear density distribution as well? 
It turns out that the answer is yes. In other words, we derive a criterion to assess the 
quality of a nuclear density distribution p{r) by observing an increase of S on going 
from crude to more sophisticated models of the nucleus. 



3 Numerical results for S 

It is well known that the kinetic energy density t is not uniquely defined. Thus, there 
are various expressions for this quantity ||13|| : 



For simplicity we have chosen the Thomas- Fermi- Weizsacker kinetic energy expres- 
sion fl3f : 

tTFW = tTF + ^X- — (12) 

72 p 

This choice allows us to use various functional forms of the density distributions, where 
no knowledge of the contribution to the density from each orbit is necessary. This is 
also corroborated indirectly by the fact that in |13|], using (12), there is not a single 
exception to the trend of the information entropy to increase as the quality of the 
electron wavefunction increases. 

We calculate S employing as an input several nuclear (matter) density distributions: 
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1) Uniform (UN) specified by the well-known rule: R = r$A 




P(r) = r° lZ o (13) 



2) Trapezoidal (TR) gfj: 

p r < c — z 

P( r ) = { Po( c + z ~ r)/2z c—z<r<c+z (14) 

r > c + z 

where c is the half-way radius and 2z is the width of the surface region. 

3) Two parameter Fermi (FM) ||17|| : 



= n (15) 

1 + exp [^) 

For the half density radius R, instead of the well known expression R = tqA" 1 ^ 3 an 
expression with an A-dependent radius parameter r = r (A) is used, which results 
from the solution of an algerbraic third order equation for the normalization of the 
density fl7|l: 



r (A) 



2 l/3 



^0 < 



1 + 



1 + 



4 / Tea 



27 Vr AV3 



1/2- 



1/3 



( 



ira 



27 \r A 1 / 3 



1/2' 



1/3' 



(16) 



where 

r = (3/4vrp ) 1/3 (17) 

It is noted that the same expression for the half density radius has also been used 
for the trapezoidal distribution (c = R) , where na = z. In both cases the parameters 
were determined by a global least squares procedure to the available experimental rms 
radii. 

4) Harmonic Oscillator (HO) Shell Model, specified by the rule : hto = 41A -1 / 3 



5) Semiphenomenological single-particle density (SP-D) of Gambhir and Patil (18 
We use the expression of ref 
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Pi 



(r-R) ~(r + R) 

l + A[l + (7^-) 2 Ne^+e^H 
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where 

A = + 

(in order to idendify R with the half-way radius), the index i = n (for neutrons) or p 
(for protons), R is a measure of the size of the nucleus and tij and oij are given in terms 
of the separation energy e« of the last particle (neutron or proton) through: 

" (19) 



2\/2mej ' 

q I m , , 

^ = ^/^ + l (20) 



Here m is the nucleon mass and q = for neutrons and g = Z — 1 for protons. 

The above Pi(r) correctly incorporates the following two important physical require- 
ments: 

a) The small r (r —>■ 0) behaviour which implies that the density contains only even 
powers of r. 

b) The asymptotic behaviour 

Pt(r) — > r~ 2ai exp(—r / di) (r — > oo) (21) 

In the present work we employ the matter distribution p(r) = p n (r) + p p (r), where 
p n (r) is normalized to the number of neutrons and p p (r) to the number of protons , so 
that p(r) is normalized to the number of particles A (mass number). 

6) Densities calculated within the framework of microscopic mean field approaches. 
Namely, densities obtained from Hartree-Fock calculations with density dependent 
forces of Skyrme type (DDHF) JT1] and the ones using the relativistic mean field 
(RMF) theory @. 

We have computed the information entropy S using the above densities. Our results 
are shown in Table 1. For the HO potential results are provided only up to 40 Ca, since 
this model is not expected to give realistic results for heavier nuclei. The results for the 
Skyrme and RMF theories have been obtained with the effective forces SkM* p^] and 
NL3 [23[] respectively. It is interesting to note, however, that using different Skyrme or 



RMF parametrizations, we have obtained quite similar results i.e. the estimate of the 
entropy is independent of the used effective force in either theory. 

We observe in Table 1 that S increases from left to the right i.e. on going from 
crude to more sophisticated (microscopic) models, although this trend is less clear to 
the right of the table, where the value of S tends to saturation. Therefore, we can 
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conclude that the larger S, the better the quality of p(r). It is also seen that from 
the various functional forms used for the description of the densities, the Fermi density 
shows the best behaviour leading to values for S quite similar to those of the microscopic 
models. It is also observed that SP-D gives also large S values i.e. it is a good density 
according to the present criterion. This is expected because SP-D reproduces fairly 
well the experimental data. Moreover, the S values of SP-D are quite similar to those 
of Fermi distribution and one could say that they both lead to p(r) of similar quality. 
The interesting thing is that the experimental input used for the determination of their 
parameters is quite different. In Fermi distribution the parameters are treated as free 
parameters determined by a global fit to experimental rms radii while in SP-D the 
parameters (not free) are adjusted to the experimental separation energies of the last 
proton and neutron of each nucleus. 



4 An alternative approach 

Gadre and Bendale (GB) 0, |J employed a different maximum entropy procedure for 
atomic systems. They studied the entropy sum S r + Sk where 

S r = -J p(f)lnp(f)df (22) 

S k = - J n(k)lnn{k)dk (23) 

S r is the position-space information entropy and Sf. is the momentum-space information 
entropy. p(f) and n(k) are the single-particle densities in configuration and momentum 
space respectively. They found that the sum S r +Sk increases with the basis-set quality 
of the electron density p(r). 

In H we have focused on the same sum S r + Sk but for another many-fermion 
system i.e. the nucleus. We have employed two models of the nucleus, first the nuclear 
mean field was approximated by a HO potential and then the model was extended by 
introducing some sort of short-range correlations. In both approaches we have found 
that S r + Sk is independent of the single parameter of the model (the HO parameter) 
i.e. it is scale invariant and characterizes every nucleus. We have also found that 
short-range correlations increase slightly the net information content S r + Sk of the 
nucleus. The relative increase in S r + Sk due to short range correlations ranges from 
2.5% for 4 He to 1.2% for 40 Ca. 

In the present work we calculate the sum S r + <Sfe for the three doubly magic nuclei 
4 He, 16 and 40 Ca using as input p(r) and n(k) from a phenomenological method [24 
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based on the natural orbital representation to construct the one-body density matrix 
(OBDM) p(f,f'). According to this method: 



p{r) = p(r,r) =J2\ i \ilj i (r)\ 2 

i 

n(k) = n{k,k) = Y J MUk)\ 2 



(24) 



where % = nlm, ipi{k) is the Fourier transform of ipi{r) and X n i are the occupation 
numbers. The natural orbitals can be looked for in the form: 



il>i(r) = ipnim(r) = R n i{r)Yi 



lm \ 



(25) 



The radial part of the natural orbitals of the OBDM are expanded |24| in terms of 
single-particle wavefunctions {4> n i{f)} preserving all the usual symmetries for spherical 
nuclei: 



(26) 



In ref. [24| three sets of s.p. wavefunctions {4> n i( r )} were used corresponding to: 

1) the harmonic oscillator potential (HO) 

V(r) = -V + ^mu 2 r 2 (V > 0) 

2) the square well (SW) potential with infite walls 

v t r ) = I ~ V ° r<x, V o >0 
I oo r > x 



(27) 



(28) 



(29) 



The parameters of the above expressions were obtained in p4| by fitting simultaneously 
the local density p(r) and the momentum distribution n(k) to the available empirical 
data. However, while for the density distribution a lot of information is available, 
this is not the case with the momentum distribution where the information is rather 
limited. Experimental data for n(k) is available only for 4 He [26]. For this reason, the 
theoretical estimates of n(k) in 16 and 40 Ca obtained within the Jastrow correlation 
method (JCM) [^7| were used in the fit. In all cases examined the best-fit values were 
obtained using the square-well single-particle wavefunctions. 



3) the modified harmonic-oscillator (MHO) potential |25 

1 B 

V(r) = -V + -mu 2 r 2 + — , V > 0, B > 

2 r z 
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5 Numerical results for the sum S r + Sk 



In table 2, we list for comparison the S r , Sk and S r + Sk values calculated by inserting 
(24) into equations (22), (23) for each set of basis wavefunctions (HO, SW, MHO). It is 
observed that S r + Sk for 4 He is a maximum for the SW case, where also the best fit to 
the experiment is obtained. However, this is not the case for 16 and 40 Ca nuclei where 
S r + Sk becomes maximum for the HO case. It should be noted, however, that for 16 
and 40 Ca theoretical input for n(k) has been used in the fitting procedure. Moreover 
the theoretical n(k) values obtained within JCM were calculated using as only input 
experimental knowledge of the density distributions. 

Our results suggest that when experimental data ( 4 He) are reproduced well, the 
sum S r + Sk is indeed a maximum, i.e. the larger S r + Sk, the better the quality of 
p(r). This also indicates that more experimental information for n(k) is called for. It 
is also seen that in 4 He where n(k) is deduced from experiment the Sk value is larger in 
SW case than in the HO and MHO cases where the values are smaller and very close 
to each other. 

6 Discussion and conclusions 

We have presented two methods of assessing the quality of nuclear density distributions. 
The first method is based on the information entropy S associated with phase-space 
distribution. It is seen that the Fermi distributions and the densities deduced from 
microscopic models correspond to larger information entropies i.e. they are of better 
quality compared to more crude models of the nucleus. It is also observed that the 
semiphenomenological density (SP-D) derived by reproducing the separation energies 
of the last nucleon and the correct asymptotic behaviour of p(r) is of comparable 
quality (according to the present criterion) with Skyrme DDHF and RMF approaches. 
This is expected because SP-D reproduces well the empirical values for the rms radii. 
Another observation is that different parametrizations of Skyrme and RMF models give 
similar results for the information entropy S. It seems that after a certain degree of 
sophistication of a nuclear model is reached, the information entropy no longer increases 
but saturates. However, such a conclusion needs more examination. 

The second method examines the sum S r + Sk of the entropies in position and 
momentum space. We have used a certain phenomenological model which constructs 
the one-body density matrix p(r,r ). The parameters of the corresponding densities 
p(r) and n(k) were found by fitting simultaneously the experimental data of p(r) and 
n(k). For n(k), however, only for 4 He some experimental input is available. Therefore, 
the predictions of a model (JCM) were used as an input for the rest of the cases ( 16 
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and 40 Ca). Our study shows that for 4 He and for wavefunctions corresponding to the 
best quality fit (SW), among those considered, the sum S r + Sk gets its highest value. 
Moreover, in 4 He, where experimental input for n(k) was used, the value for the best 
quality wavefunctions is significantly larger compared with the values obtained with the 
use of the other wavefunctions. For 16 and 40 Ca , however, the sum becomes maximum 
for HO wavefunctions which do not correspond to the best fit. This should be connected 
with the use of theoretical predictions for n(k) due to lack of experimental information. 
These results, however, should be considered preliminary because a limited number of 
cases have been examined in the OBDM method of ref. || , though they are in the 
correct direction and suggest that the MEP works, in accord with the first method of 
the present paper. 



Table 1: The information entropy S for several nuclei using various nuclear distributions 
(GBP approach) 



Nucleus 


UN 


TR 


HO 


SP-D 


FM 


DDHF 


RMF 


16 o 


99.69 


101.68 


102.06 


102.11 


102.12 


102.10 


102.13 


32 S 


199.38 


202.34 


202.75 


202.81 


202.88 


202.80 


202.90 


40 Ca 


249.23 


252.59 


253.04 


253.08 


253.18 


253.11 


253.17 


90 Zr 


560.77 


566.10 




566.66 


566.96 


566.90 


566.94 


116 Sn 


722.76 


728.68 




729.45 


729.91 


729.97 


729.88 


208pj 3 


1295.99 


1304.57 




1304.91 


1305.99 


1305.81 


1306.20 



Table 2: The information entropies <S r , Sk and S r + Sk for three nuclei using three sets of 
single particle wavefunctions corresponding to harmonic oscillator (HO), square well(SW) 
and modified harmonic oscillator (MHO) potentials (GB approach) 



Nucleus 




HO 






SW 






MHO 






<S r 


s k 


S r + Sk 




Sk 


S r + Sk 




Sk 


S r + Sk 


4 He 


7.27 


2.06 


9.33 


6.91 


4.26 


11.17 


7.10 


2.28 


9.38 


16 o 


27.96 


3.69 


31.65 


27.91 


2.80 


30.71 


28.02 


1.20 


29.22 


40 Ca 


66.22 


7.82 


74.04 


65.67 


-5.81 


59.86 


66.38 


-9.09 


57.29 
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